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ABSTRACT 

Aims. The discovery of a statistically significant clustering in the distribution of gamma-ray bursts (GRBs) has recently been reported. 
Given that the cluster has a characteristic size of 2,000-3,000 Mpc and a redshift between 1.6 < z < 2.1, it has been claimed that 
this structure is incompatible with the cosmological principle of homogeneity and isotropy of our Universe. In this paper, we study 
the homogeneity of the GRB distribution using a subsample of the Greiner GRB catalogue, which contains 314 objects with redshift 
0 < z < 2.5 (244 of them discovered by the Swift GRB Mission). We try to reconcile the dilemma between the new observations and 
the current theory of structure formation and growth. 

Methods. To test the results against the possible biases in redshift determination and the incompleteness of the Greiner sample, we 
also apply our analysis to the 244 GRBs discovered by Swift and the subsample presented by the Swift Gamma-Ray Burst Host 
Galaxy Legacy Survey (SHOALS). The real space two-point correlation function (2PCF) of GRBs, ^(r), is calculated using a Landy- 
Szalay estimator. We perform a standard least-;^^^ fit to the measured 2PCFs of GRBs. We use the best-fit 2PCF to deduce a recently 
defined homogeneity scale. The homogeneity scale, R/j, is defined as the comoving radius of the sphere inside which the number of 
GRBs N(< r) is proportional to r’ within 1%, or equivalently above which the correlation dimension of the sample D 2 is within 1% 
ofD 2 = 3. 

Results. For a flat ACDM Universe, a best-fit power law, = (r/ro)^’', with the correlation length ro = 413.64 ± 135.40 ti^'Mpc 
and slope 7 = 1.57 ± 0.63 (at Icr confidence level) for the real-space correlation function f(r) is obtained. We obtain a homogeneous 
distribution of GRBs with correlation dimension above D 2 = 2.97 on scales of r > 8,200 /r'Mpc. For the Swift subsample of 
244 GRBs, the correlation length and slope are ro = 387.51 ± 132.75 /j^'Mpc and 7 = 1.57 ± 0.65 (at Icr confidence level). The 
corresponding scale for a homogeneous distribution of GRBs is r > 7,700 /r^'Mpc. For the 75 SHOALS GRBs, the results are are 
ro = 288.13 ± 192.85 /r^'Mpc and 7 = 1.27 ± 0.54 (at Icr confidence level), with the homogeneity scale r > 8,300 /r^'Mpc. For the 
113 SHOALS GRBs at 0 < z < 6.3, the results are ro = 489.66 ± 260.90 /r'Mpc and 7 = 1.67 ± 1.07 (at Icr confidence level), with 
the homogeneity scale r > 8,700 /i^'Mpc. 

Conclusions. The results help to alleviate the tension between the new discovery of the excess clustering of GRBs and the cosmo¬ 
logical principle of large-scale homogeneity. It implies that very massive structures in the relatively local Universe do not necessarily 
violate the cosmological principle and could conceivably be present. 
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1. Introduction 


Gamma-ray bursts 
in our Universe, 
studied by 


(GRBs) are the most 
Their c o smolo gi cal 

(I 1973 I): iMeegan et al.l (Il992l): 


energetic events 
origin has been 


Klebesadel et al 


Kouveliotou et al.l (Il993l): ICostaet^ (Il99^: iParadi 


il997l) : iHarrison et al.l (Il999l) : iMeszaros & Gehrelsl ( 2012I) . 
GRBs and luminous red galaxies (LRGs) are both luminous 
tracers of matter in our Universe. Unlike most LRGs, GRBs 
have a larger redshift range that reaches up to z ~ 8. They 
have minimum separations over 100 Mpc. Therefore, they 
are valid indicators of potential large-scale structures in the 
intermediate-redshift Universe. 

According to modern cosmology, structures on different 
length scales in our Universe all have their origins in the quan¬ 
tum fluctuations of the inflation field, the scalar held which gen¬ 
erates the inflation after the birth of our Universe. These pri¬ 
mordial Gaussian random phase fluctuations, which later lead 
to the density fluctuations in different modes, enter the horizon 


on different epochs and grow as time passes by. This gives rise 
to the hierarchical scenario of the matter clustering - the fluc¬ 
tuations with longer comoving wavelengths enter the horizon 
later and have less time to evolve than the small-scale fluctua¬ 
tions. As large-scale structures are closely related to these long- 
wavelength modes, one can conclude that a finite age of our Uni¬ 
verse would result in a limited maximum size of the large-scale 
structure in our Universe. Translated into the language of modern 
cosmology, it is the cosmological principle - the matter distribu¬ 
tion in our Universe is homogeneous and isotropic on sufficiently 
large scales. On small scales, our Universe is inhomogeneous, 
with structures like galaxies and galaxies clusters. 

The transition scale between homogeneity and inhomogene¬ 
ity is called the ‘homogeneity scale lI The homogeneity scale 
of the distribution of ma tter has lon g been studied and the re¬ 
sults are quite scattered. iHogg et ^ (l2005h analysed the enor¬ 
mous LRG samp le of the Sloan Digital Sky Survey (SDSS) 
(lYork et al.ll200C)l) and presented a homogeneity scale Rh ~ 70 


* limh@ihep.ac.cn 
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’ Alternatively, one can call it the in homogeneity scale. Here we use 
the conventional term that was used in lYadav et al.l (l2010h . 
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/i~* Mpc. Similar results we re obtained by ISarkar et alJ (l2009h 
and lScrimgeour et"^ (l2012h . who performed a multifractal anal- 
ysis over 200,000 blue l u minous galax i es in t he WiggleZ survey 
ilDrinkwater et al.l|2010h . iLabini et al.l (l2009h claimed to find a 
homogeneity scale above 100 /z"'Mpc after studying the galaxy 
cat alogue of SPS S. 

iHorvath et al.l (l2014t) have recently reported the discovery 
of a statistically significant clustering in the GRB sample be¬ 
tween 1.6 < z < 2.1. They called it the Hercules-Corona Bo¬ 
realis Great Wall (Her-CrB GW). It has a characteristic scale 
of ~ 2,000 Mpc and its longest dimension ~ 3,000 Mpc, is 
six times larger than the size of the Sloan Great Wall. Its char¬ 
acteristic size is far above the upper limit of the homogeneity 
scale placed by the fractal dimensional analysis based on the 
galaxy surveys. Th e two-dimensional Kolmogorov-Smirnov test 
(iLones et al.ll2008h shows a 3cr deviation. The clustering excess 
cannot be entirely attributed to the known sampling biases. The 
existence of the potential structure, defined by the GRBs, is con¬ 
sidered to be inconsistent with the cosmological principle and 
beyond the standard excursion set theory for structure growth. 

In this paper, we examine the spatial distribution of the GRBs 
using a subsample of Greiner’s GRB catalogue, which contains 
314 objects with redshift 0 < z < 2.5 (244 of them discovered 
by the Swift GRB Mission). The sample encompasses the red- 
shift region of the reported potential structur e mapped by the 
GRBs. We first use the Landy-Szalay estimator (iLandv & Szalavl 
Il993h to estimate the real space two-point correlation function 
(2PCF) ^(r) of the GRB sample. We fit a simple power law to 
the measured GRB ^(r). We then use the best-fit f(r) to deduce 
the correlation dimension £> 2(0 and the homogeneity scale Rh 
of the GRB distribution. To test the results against the possi¬ 
ble biases in redshift determination and the incompleteness of 
the Greiner sample, we also apply the analysis to the subsample 
presented by the Swi ft Gamma-Ray Burs t Host Galaxy Legacy 
Survey (SHOALS) in iPerlev et al. I (l2015h . The results are plot¬ 
ted in the same figure for comparison. 

The rest of the paper is organized as follows. In Section 2, 
we introduce the GRB sample and the techniques used in our 
analysis. The real-space 2PCF ^(r) of the GRB sample and its 
best-fit power law are both obtained in this section. In Section 
3, the definitions and calculations of homogeneity scale and cor¬ 
relation dimension are presented. We deduce the specific Rh for 
the GRBs with 0 < z < 2.5. Relevant physical implications and 
comparisons to the results of other surveys are presented in Sec¬ 
tion 4. 


2. Data and techniques 

2.1. GRB catalogue and the subsample 

In this work we primarily use a sample of 314 GRBs with red¬ 
shift 0 < z < 2.5. All of these GRBs are from the collec¬ 
tion presented by J. Greiner (2014j3; 244 of them come from 
the NASA Swift Mission; and the rest come from BeppoSAX 
GRBM, HETE2, IPN, and INTEGRAL, etc. We use the data re¬ 
leased on July 8, 2015. The entire sample contains more than 
1,000 objects. Only 431 of them have well-measured redshifts. 
The well-measured subsample has a redshift 0 < z < 9.2. In the 
redshift region z > 7, there are only two GRBs: GRB 090423 
(z = 8.26) and GRB 090429B (z = 9.2). They are omitted from 
our numerical analysis since they are of little statistical signif¬ 
icance. This leaves a subsample of 429 GRBs at 0 < z < 6.7. 


^ http://WWW.mpe.mpg.de/~j cg/grbgen.html 


Given that the potential GRB structure reported bv iHorvath et al.l 
(l2014l) has a redshift 1.6 < z < 2.1, we cut this sample down fur¬ 
ther by using only those GRBs having well-determined redshifts 
at 0 < z < 2.5, which encompasses the redshift region of the 
potential structure. We have a final subsample of 314 GRBs at 
redshift 0 < z < 2.5. 

The redshifts of GRBs can be estimated in a number of 
ways, i.e. through the absorption spectroscopy of the optical af¬ 
terglow (the vast majority of the Swift sample are measured in 
this way) or measuring the emission lines of their host galaxies 
(which is observationally expensive). Eaintness of the afterglows 
of some events have left the sample intrinsically optically biased. 
Alter native approaches have been proposed to solve this prob- 
lossi I ' " ' 


iHunt et al. 2014). One of them is to develon a series of observ- 
ability cuts with a set of optimized parameters to isolate a subset 

of the GRB sample (Jakobsson et al. 

2006t ICenko et alJl2006l 

iPerley et al. 11200^ iGreiner et al. 2011 

). Using this method, one 

can obtain a GRB sample whose afterglow redshift completeness 

is close to about 90% (Hiorth et al.ll2012; 

Jakobsson et alJ 20121 

iKruhler et al.ll2012t ISalyaterra et al. 2012 

; Schulze et al.ll2015l). 


the level that is necessar y for systematic b iases not to domi¬ 
nate t he statistical ones JPerley et al. ll20T5l) . Using this tech¬ 
nique, |Periex£tar] (1201 5h established the largest and most com¬ 
plete (92% completeness) GRB redshift sample to date, which is 
called the SHOALS sample. The SHOALS sample contains 119 
objects in total, 112 of which haye well-determined redshifts at 
0 < z < 6.3 (75 are at 0 < z < 2.5). To test the results against 
the possible biases in redshift determination and the incomplete¬ 
ness of the Greiner sample, we also apply our analysis to this 
subsample of GRBs. 

Thus, we base our analysis primarily on the Greiner sam¬ 
ple of 314 GRBs. In comparison, we also apply our analysis 
to the Swift subsample of 244 GRBs and the SHOALS yery 
complete GRB sample. The angular and redshift distributions of 
these samples are plotted in Eigure [Ta] and [10 The correspond¬ 
ing celestial coordinates and the redshifts of the objects in the 
Greiner sample and the SHOALS sample are listed respectiyely 
in Tables 1 and 2, which are publicly ayailable onlin^l Most 
of them haye a comoying separation r > 100 /r 'Mpc, with the 
longest separation distance up to r ~ 10 /r *Gpc. 


2.2. Two-poInt correlation function 


Giyen a GRB in the spatial yolume dVi, the 2PCE of GRBs, f(r), 
is defined as the probability of finding ano ther GRB in dV 2 with 
a separation distance r, i.e. (lPeebleslll980h 

dPn = «^[1 + ^(r)]dVidV2, (1) 


where n is the mean number density of the GRBs. To calculate 
^(r), an auxiliary random sample of Nr points is generated in a 
window W of obseryations. A window VT is a three-dimensional 
space of yolume V, the same yolume as that on which the obser- 
yation was made. 

A statistical estimation of ^(r) inyolyes a pair count of neigh¬ 
bouring GRBs at a giyen separation scale. The most widely 
used estimator of the 2P CP is the Landy-Szalay estimator |ls(^) 
(iLandy & Szalay|[T993h . 


|ls(G = 


DD{r) - 2DR{r) + RR(r) 
RR{r) 


( 2 ) 


^ Tables 1 and 2 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 
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90°N 



(a) Anguair Distribution 



(b) Redshift Distribution 


Fig. 1. Angular and redshift distributions of the GRB samples, (a) The 
angular distribution of the GRB samples in J.2000 equatorial coordi¬ 
nates. The red solid dots represent the 244 GRBs at 0 < z < 2.5 detected 
by Swift, while the black solid dots represent those discovered by other 
detectors within the same redshift range. The red and black solid dots 
constitute Greiner’s GRB sample of 314 objects at 0 < z < 2.5. The 
blue circles represent the 112 GRBs (at 0 < z < 6.3) from SHOALS, 
(b) The redshift distribution of the GRB data. The i/-axis denotes the 
number of objects in each redshift bin. The green shaded area plus the 
purple area indicates the total of 314 GRBs from the Greiner sample. 
The dashed line represents the contribution from the SHOALS subsam¬ 
ple of 75 objects at 0 < z < 2.5. 


where DD{r) and RR{r) are, respectively, the number of GRB 
pairs within the seperation d e{r — hrjl, r + ts.rjl} (Ar is the bin 
width used in the statistical estimation of f(r)) in the observed 
data set D and in the auxiliary random sample R in the window 
W, while DR(r) is the number of GRB pairs between the ob¬ 
served data and the random sample with the same separation. 
The parameter d is the comoving seperation distance of GRBs. 
Specifically, DR(r) is defined as DR{r) = NDR{r)/{NDNR) {No 
and Nr are, respectively, the total number of GRBs in the data 
set D and in the rand om sample R), where NoRir) is given by 
(iKerscher et alJl2()()(ih 


NoRir) - 

xeD yeR 


(3) 



separation h ' Gpc 


Fig. 2. The distribution of the comoving separations d between the 
GRBs in the Greiner sample (including the Swift subsample of 244 
GRBs and most of the SHOALS GRBs). The x-axis is the comoving 
separation of the GRB data, in units of /t^'Gpc. The ^-axis is the pro¬ 
portion of the GRB number in each bin to the whole sample, which has 
been normalized to 1. The separation d obeys a Gaussian distribution 
with the expectation value J 3,710 /t^'Mpc. 


The summation runs over all the coordinates of GRBs (repre¬ 
sented by X and y) in the observed data set D and the random 
sample R in the window W. The value of the function F{x, y) 
equals 1 when the separation of the two objects is within the 
distance d(x, y) e [r - Ar/2, r + Ar/2] or otherwise equals 0. 
The expressions DD{r) = A^z)D(r)/[A^D(A^D - 1)] and RR{r) = 
NRR{r)l{NR{NR - 1)] respectively represent the normalized num¬ 
ber of GRB pairs within the separation mentioned above in the 
observed data set D and in the random sample R; Nooir) and 
NRR{r) are defined in a similar way as NoRir) in equation Q. 

We use the jackknife resampling method to determine the 
statistical uncertainty of the measured 2PCF of GRBs. The 
jackknife method is an internal method of error estimation 
that is extensively us ed to determine the errors of 2PCF of 
galaxies and quasars (iRoss et al.l 120071: ISawangwit et al.ll201 ll: 
iNikoloudakis et al.ll2013l) . The entire sample is divided into N' 
subsamples of roughly equal size. The jackknife error estimator 
is given as 


N' 


o-Lk('-) = Z 
(■'=1 


DRi'jr) 

DR(r) 




(4) 


where (r) denotes the estimate of the 2PCF on all of the {N' -1) 
subsamples except the /-th one. 


2.3. Calculating the two-point correlation function, ^(r) 

The estimation of the 2PCF of GRBs, f(r), is made by counting 
the pairs in and between the observed GRB distribution and a 
catalogue of randomly distributed GRBs. The density of random 
points that we use for the estimation is 20 times the density of 
the GRB data. We use a bin width of Alog[g(r) = 0.2 'Mpc. 
The estimation of f(r) depends on the assumed cosmology. We 
use a flat ACDM cosmological model (cold dark matter plus a 
cosmological constant A) in this work, with Qa = 0.72, Qm = 
0.28, Ho = 100/2 km s^' Mpc^*, h = 0.7. 

Furthermore, the 2PCF is usually measured in redshift space. 
As shown in Figure|2] the GRBs in the observed sample have an 
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Fig. 3. Best-fit power law of the measured real-space 2PCF ^(r) at 
200h^' < r < 10‘'/j“' Mpc. We assume a flat ACDM cosmological 
model, with £2 a = 0.72, = 0.28, Hq = 100/r km s“*Mpc“', h = 0.7. 

The real-space 2PCF ^(r) measured via equation (El for the Greiner, 
Swift, and SHOALS samples are, respectively, indicated by green cir¬ 
cles, red solid triangles, black solid diamonds, and cyan squares with 
Icr jackknife error bars that are estimated from 0. The density of 
random points we use for the estimation is 20 times the density of 
GRB data, ^(r) is plotted in equally spaced logarithmic intervals of 
Alogjofr) = 0.2 /i^'Mpc. The best-fit power laws of the form 0 for 
the measured f(r) are plotted in solid lines, with the best-fit parameters 
given in the legends. 


average separation distance over 100 /t^'Mpc (with the largest 
separation up to ~ 10 /t 'Gpc). On such a scale, the current struc¬ 
ture formation theory predicts that the evolution and clustering 
of matter should still be in the linear regim e at the present time 
(ISoringel ^ al.ll2()()5t lEisenstein et al.ll2007ll . The redshift-space 
distortion^ due to the small-scale peculiar velocities of the ob- 
j ects and the reds hift variances are also minimal on this scale 
(IRoss et al.l l2007h . Thus, the difference between the redshift- 
space and the real-space correlation functions on such large 
scales could be negligible. For convenience, the calculation and 
analysis in this paper are done in real space. 

For the estimation of the jackknife error we take 

N' - 5, and split the sample into five redshift regions with 
equal redshift intervals Az - 0.5. The real-space 2PCFs in a flat 
ACDM Universe for Greiner’s 314 GRBs, the subsample of 244 
objects discovered by Swift, and the SHOALS subsample are re¬ 
spectively plotted in Figure[3]in equally spaced logarithmic inter¬ 
vals. It spans a comoving distance of scale 10^/i“' < r < 

Mpc . 

A power law of the form 

(5, 


is usually fitted to the correlation functions of g alaxies 
and galaxy cluster s (|Davis_&Peeb^ 19831: Bahcalll 
Maddoxet^ 19901: iPeacock & West Il992t iDaltonet akl 


1988; 


1994; 


Zehavi et al.l 120041) . The parameter ro is the comoving correla¬ 
tion length, in units of /r *Mpc. The slope y is a dimensionless 
constant. Similarly, we fit a power law of the form in Eq. (|5]l to 
the measured 2PCF data over the range 200/2 ' < r < 10^/2 ' 


A more fo rmal descr i ption of ‘redshift distortions’ can be found in 
Section 9.4 in lDodels^ d2008h . 


Mpc. We perform a standard least-;^^ fit. For a flat ACDM Uni¬ 
verse, the best-fit values of the correlation length ro and slope 
y for the real-space correlation function f(r) for Greiner’s 314 
GRB data at z < 2.5 are (with Icr confidence level errors) 


ro = (413.64 + 135.40) /i^'Mpc, (6) 

y= 1.57+ 0.63, (7) 

with;^^^j^ = 0.73f. For the subsample of 244 objects discovered 
by Swift, we obtain (with Icr confidence level errors) 

ro = (387.51 + 132.75) /i^^Mpc, (8) 

y= 1.57+ 0.65, (9) 

withuL = 0.90. 

For the SHOALS subsample of 75 GRBs, the results are 
ro = (288.13+ 192.85)/r^'Mpc, (10) 


y= 1.27+ 0.54, (11) 

W\\h - 0.50. For the SHOALS sample of 113 GRBs at 
0 < z < 6.3, the results are 

ro = (489.66 + 260.90) /z^'Mpc, (12) 

y= 1.67+ 1.07, (13) 

with . = 2.82. From these results, one can see that the best-fit 
values of ro and y increase with the growing number of GRB data 
points as well as higher redshifts. We plot the best-fit power-law 
models of ^(r) in Figure [3] 


3. Homogeneity of the GRB distribution 

Given the best-fit 2PCF, ^(r), for the GRB sample in the previous 
section, we are now in a position to calculate the homogeneity 
scale, Rh, for the GRB distribution at 0 < z < 2.5. We first give 
a brief introduction of the correlation dimension D 2 {r) for a ran¬ 
dom distribution of data points. We describe the relation between 
the value of D 2 (r) and the concept of a homogeneous distribu- 
tion. A more form al tr eatment of this section can be found in 
lYadav et'aP (1201 Oh and lScrimgeour et"^ (l2012h . 

3.1. Correlation dimension, D 2 (r) 

Several methods have been developed to investigate the homo¬ 
geneity of the gala xy distribution. The most popular among them 
is fractal analysis dYadav et al.ll2005|l . A fractal is a kind of ge¬ 
ometrical object where every small part of it appears as a re¬ 
duction of the entirety. In fractal analysis, the concept ‘fractal 
dimension’ is invoked to describe the homogeneity of the dis¬ 
tribution of a point set. One of the most common definitions 
of fractal dimension is the ‘correlation dimension’, D 2 (r). Un¬ 
like other homogeneity indicators, deviations caused by a size- 
limited sample would only result in second-order changes to 
D 2 {r). Thus, D 2 {r) is regarded as a robust measure of homogene¬ 
ity and is extensively used in the homogeneity investigations of 
galaxies and qu asars (iLabini et al.l2009t;IScrimgeour et al.ll2012l : 
lNadathuil2ni'l . 

In this paper, we use the working definition of D 2 {r) given in 
IScrimgeour et al.l (l2012h to study the homogeneity of the GRB 
distribution. We limit our discussions to a three-dimensional 
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space. Given a set of points in space, the measurement of D 2 {r) 
is to find the average number of neighbouring points, N{< r) 
inside a three-dimensional sphere of radius r centered at each 
point. One can formulate the scaling behaviour of N{< r) as 

N(< r) oc (14) 


where D is the fractal dimension of the distribution. A^(< r) can 
actually be given as N{< r) - AurTilh, where n is the mean 
number density of points in that region. For a homogeneous dis¬ 
tribution of the point set, h is a universal constant and N{< r) 
then scales as oc r^. From equation (fT4l) . this implies D - 3. An 
inhomogeneous distribution would result in a non-universal h. 
From equation (fT4li . this gives D < 3 or D > 3: D < 3 stands 
for an inhomogeneous distribution of the point s et, whileD > 3 
repre sents a ‘super-homogeneous’ distribution (iGabrielli et al.l 
l2002h . In the literature, N{< r) is usually divided by the number 
expected for a homogeneous distribution, Anr^hl3 (here n is a 
universal constant), to correct for incompleteness. The corrected 
N{< r) is given as 

N{< r) cc r^-^. (15) 


For a homogeneous distribution, one has 7V(< r) oc - 1. 

In general, the correlation dimension D in equation (fT4l) is a 
function of the sphere radius, r, i.e. D - D 2 {r). It can be deduced 
from the count-in-sphere number N{< r) as 


£>2(r) = 


d lnA(< r) 
d Inr 


d lnAf(< r) 
d In r 


(16) 


The correlation dimension D 2 {r) measures the scaling properties 
of N{< r) without being affected by the amplitude of N(< r) 
(which is related to the mean number density of the data points 
in the regions.) It is therefore an objective measure of the homo¬ 
geneity in the statistical analysis of matter distribution. To de¬ 
duce the correlation dimension D 2 (r) for the GRB distribution, 
one should first calculate N{< r). For the GRBs N{< r) can 
actually be obt ained by integrating the 2PCF ^(r) for the GRBs 
(iPeehleslIlQSntl : 


f [1 +^(r')]47rr'W. (17) 

47Tr^ Jo 

Combining equation (fThl l with (fTTl i. we obtain an analytical 
expression of D 2 {r): 


D2{r) - 


d Inr 


In 


Anr^ 


f [\ + ^{r')]Anr'~dr' 

Jo 


-h3. 


(18) 


Given the best-fit power law of ^(r) as in equation (|5]l with the 
correlation length ro and slope j of ^(r') given in equations (|6]l 
to (foi l, we are now in a position to estimate the homogeneity 
scale of the GRB distribution. 



r [*-' Mpc] 


Fig. 4. Correlation dimension, Diir), for different GRB samples. The 
D 2 {r) calculated from equation l ll8i is shown in solid lines. The dashed 
line indicates the critical value defined for the transition from a homo¬ 
geneous to an inhomogeneous distribution of the GRBs, i.e. 1% from 
the homogeneity, Diir) = 2.91. 


is not affected by the sample size. Given the correlation dimen¬ 
sion £> 2(0 of the GRB distribution, the homogeneity scale Rh is 
defined as the scale on which D2(r) of the sample is within 1% 
of D2 = 3, i.e. D2{r = Rh) - 2.91 . An equivalent definition of 
Rh can also be given by N{r), for which Rh is defined as the co¬ 
moving radius r of the sphere inside which the number of LQGs 
N(< r) is proportional to P within 1%, i.e. N{r - Rh) - 1.01. 

Using these definitions, we calculate the homogeneity scale 
of the distribution of the GRB sample. The theoretical predic¬ 
tion of D2{r) and the corresponding value of Rh for the GRB 
samples are plotted in Figure |4] for comparison. For Greiner’s 
sample of 314 GRBs, the homogeneity scale for the GRB dis¬ 
tribution is Rh — 8,200 h *Mpc, which means that the cos¬ 
mological principle is retained on such a scale (r > Rh). One 
should observe a homogeneous distribution of GRBs on scales 
r > Rh. For the subsample of 244 GRBs discovered by Swift, 
the result is Rh - 7,700 h 'Mpc. For the SHOALS subsample 
of 75 GRBs at 0 < z < 2.5, we have Rh = 8,300 h 'Mpc. For 
the SHOALS subsample of 113 GRBs at 0 < z < 6.3, we have 
Rh = 8,700/2-'Mpc. 

The scale of homogeneity Rh for the GRBs can also be 
considered as the upper limit of the characteristic size of any 
clustering and structures in the GRB distribution. The potential 
structure mapped by the GRBs at 1.6 < z < 2.1 reported by 
iHorvath et al.l (1201 4l) has a characteristic size of ~ 2,000 Mpc 
and is thus well within the limits we obtain. Therefore, we con¬ 
clude that the existence of such an excess clustering in the GRB 
distribution is still compatible with the cosmological principle 
and the standard theory of structure formation. 


3.2. Scale of homogeneity, Rh 

There are several ways to define the homog eneity scale Rh with 
respect to the correlation dimension D 2 (r). lYadav et"^ (1201 Ol) 
defined the homogeneity scale as the scale above which the devi¬ 
ation of the fractal dimension D 2 {r) from the ambient spatial di¬ 
mension becomes smaller than the statistical dispersion of D 2 {r) 
itself. For the limited size of the GRB sample that we use in 
our work, the statistical dispersion of D 2 ir) might be large. This 
effect of a limited sized sample would therefor e enlarge the esti¬ 
mation of Rh for our GRB sample. Following IScrimgeour et al] 
(l2012h . in our analysis we use a more robust definition of Rh that 


4. Conclusions and discussion 

Redshift surveys (iDrinkwater et al.ll2010t lAhn et al']|2014ll pro¬ 
vide about a hundred thousand galaxies that can be used for a 
homogeneity investigation of the matter distribution in our Uni¬ 
verse. Most of these galaxies have redshifts z < 1. GRBs usu¬ 
ally have a larger redshift range than the galaxies (which can 
reaches up to z ~ 8) and thus provide valid luminous indica¬ 
tors of matter distribution in the intermediate-redshift universe. 
Quasars have redshifts as high as GRBs, but their observations 
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only cover a limited sky are^l This involves specific technical 
treatment in data process when using the Landy-Szalay estima - 
tor to estimate the 2PCF for quasars (iKaragiannis et al.ll2014l) . 
Compared with galaxies and quasars, the GRB sample has a 
full-sky angular distribution and a redshift range up to z ~ 8. 
In fact, in our study we found that most GRBs have a comoving 
separation > 100 h 'Mpc (see Figure|2]i and therefore they pro¬ 
vided an efficient probe of matter correlation on large scales in 
the intermediate-redshift universe. By ‘efficient’ here, we mean 
that the size of the data we used to estimate ^(r) on such a scale 
is much less than that used by galaxy or quasar surveys, since the 
number density of galaxies and quasars observed is much higher 
than GRBs. Using galaxies or quasars to probe the correlation 
and clustering on such a scale (e.g. r ~ 200 h^^Mpc) would in¬ 
volve counting many more objects inside the sphere of radius 
r and therefore would cost more computing machine time than 
using GRBs. 

In this paper we used a sample of 314 GRBs from the col¬ 
lection presented by J. Greiner (2014) to study the homogene¬ 
ity of matter distribution on large scales. They cover a redshift 
range 0 < z < 2.5, 244 of which discovered by the Swift GRB 
Mission. We calculated the real-space 2PCF, ^(r), for GRBs 
on scales lO^/r ' < z < lO^^/r ' Mpc. The measured ^(r) for 
GRBs on scales r > 200 h *Mpc can be well ht by a power 
law of the form ^(r) = (rlroy^, with correlation length ro = 
(413.64 + 135.40) h^'Mpc and slope y = 1.57 + 0.63 (Icr confi¬ 
dence level), - 0.73. For the subsample of 244 objects 

discovered by Swift, we obtained ro = (387.51 + 132.75)/2“'Mpc 
and y = 1.57 + 0.65 (Icr conhdence level), - 0.90. 

We then used the best-ht f(r) to deduce the homogeneity 
scale Rh for the GRB distribution at 0 < z < 2.5. We ob¬ 
tained Rh - 8,200 h^'Mpc. For the subsample of 244 objects 
discovered by Swift, we obtained Rh - 7,700 h^^Mpc, which 
means that above such a scale GRBs can be considered to have 
a homogeneous distribution. On scales r < Rh, an inhomoge¬ 
neous distribution of GRBs is assumed for the standard excur¬ 
sion set theory of str ucture growth. The po tential GRB structure 
recently reported by iHorvath et al.l (l2014l) has a characteristic 
size of ~ 2,000 Mpc with its longest dimension ~ 3,000 Mpc. 
Both are well below the homogeneity scale Rh we deduced for 
the GRB distribution. This implies that the discovery of such an 
angular excess of the GRB distribution is still compatible with 
the cosmological principle, which assumes that the matter distri¬ 
bution of our Universe is homogeneous and isotropic over a large 
smoothing scale. For the distribution of GRBs, we suggested 
that such a scale is Rh - 8,200 h *Mpc (which corresponds 
to Rh - 11.7 (ffic given the choice of h - 0.7). The comoving 
cosmic horizorQ within which causality holds are /h - 10 Gpc 
at z - 2 (the redshift of the potential GRB structure). The Rh we 
obtained for the GRB distribution is slightly larger than the Ih at 
z - 2. We hope that with the growing size of the observed GRB 
sample and a more accurate measurement of the 2PCF ^(r) of 
GRBs, this difference will be eliminated by future observations 
and investigations. 

To test the results against the possible biases in redshift de¬ 
termination and the incompleteness of the Greiner sample, we 
also applied our analysis to the GRB subsample presented by 
SHOALS in iPerlev et al. I (l2015h . Despite its limited sample size 
(with a total of 119 objects), the afterglow redshift completeness 


^ The quasar catalogue from the SDSS DRIO contains 166,583 quasars 
detected over 6,373 deg^. Over half of them have redshifts z > 2.15. 

® The comoving cosmic horizon is given by G = C da' where 

H(a) is the Hubble factor and a is the cosmic scale factor. 


of the SHOALS GRB sample is 92%. It was proposed in or¬ 
der to strike a balance between redshift completeness and over¬ 
all statistical size; the SHOALS GRB sample provides largest 
and most complete GRB redshift sample for clustering analy¬ 
sis to date. It has been employed to test the viability of our 
results against the biases and incompleteness of the Swift and 
Greiner samples. We obtained ro = 288.13 + 192.85 h 'Mpc and 
y = 1.27 ± 0.54 (at Icr C.L.), with;tL = for the SHOALS 
75 GRBs at 0 < z < 2.5. The deduced homogeneity scale is 
Rh = 8,300 h 'Mpc. For the 113 objects in the SHOALS sam¬ 
ple that have well-determined redshifts at 0 < z < 6.3, we ob¬ 
tained ro = 489.66 ± 260.90 h^'Mpc and j - 1.67 + L07(at 
Icr C.L.), with x^[n - 2.82. The deduced homogeneity scale 


is Rh - 8,700 /i^'Mpc. Two comments are necessary. First, the 
best-fit values of ro and y increase with growing number of GRB 
data points as well as higher redshifts. This manifests itself in 
the results obtained from SHOALS’s GRB sample, and also in 
the results obtained from the Greiner sample. Second, one can 
see that the Rh obtained from the SHOALS very complete GRB 
sample do not differ much from those obtained from the Greiner 
and Swift subsets, which implies that the Greiner GRB sample 
and the Swift GRB sample are both robust in the statistical anal¬ 
ysis of clustering and matter distribution. 

It should also be noted that for distributions of different ob¬ 
jects, the scale of homogeneity Rh are different because the dif¬ 
ferent objects have different origins and evolution dynamics, and 
thus have different bias factors. The 2PCF can also be given by 
^(r) s h^^inass(r), where b is the bias factor for specihc kinds of 
objects. Given the same 2PCF ^mass(0 of the underlying (dark) 
mattefl the amplitude of ^(r) depends on the value of b. In gen¬ 
eral, h is a function of comoving radius r and redshift z, i.e. 
b = b{r,z). Thus, for different kinds of objects with different 
characteristic (mass) scales and redshifts, b{r,z) take different 
values, which would result in a varying amplitude of ^(r). From 
equation (fTSl) . one can expect a different homogeneity scale Rh 
for the distribution of differ ent objects. For galaxies, Rh takes a 
value of ~ 70-100 /r^ 'Mpc (lYadav et al.ll2005l:lHogg et al.l|2005l: 

li. 


ISarkar et ^ l2009t iLabinil 1201 It Scrimgeour et all 


quasars in the DR7QSO c atalogue dSchneider et al 


value is Rh ~ ISO h 'Mpc (lNadathud2013h . For GRBs with red- 


u32 


20121). For 


Ma), the 


shift 0 < z < 2.5, we obtained a Rh - 8,600 h 'Mpc. Our study 
provides a supplement to t he clustering analysis based on the 
galaxie s jYadav et al.l2005l; Hogg et al.l2005t|Sarkar et al.l2009l; 


Labini 201 ll: IScrimgeour et al.(l2012l) and qua sars d White et aTl 

2012l:lNadathudl2013HKaragiannis et al.ll2014l) . 


Homogeneity on large scales is one of the cornerstones of 
modern cosmological theory. The large-scale homogeneity in 
the very early Universe (at redshift z - 1100) is well sup¬ 
ported by the high degree of isotropy of the cosmic microw ave 
background radiation power spectrum (iBennett et alj|201^ and 
the PLANCK data (lAde et al.l l2013h . Although recently the 
discoveries of several large quasar groups (LQGs) and struc¬ 
tures have been reported , e.g. the CCL QG (i.e. UL28) an d 
Ul.ll Jg^es et al.l2012h.the Huge-LOGOowes et alJ2013h . 


and the Her-CrB GW ( Horvath et al.ll2014l) . the method used 
for assessing the statistical si gnificance and ov erdensity of the 
groups varies with individuals iNadathud (l20 1 3l) provided a frac¬ 
tal dimension analysis of the DR7 quasar catalogue and found 
that the quasar distribution is homog eneous above the scale 
~ 130 'Mpc. IChangbom et al.l (120121) carried out a large cos- 


^ The underlying matter 2PCF, ^mass(r) is given by the Fourier trans¬ 
form of the primordial matter power spectrum, P(k) which is generated 
by the cosmic inflation dScrimgeour et alJ2012h . 
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mological A^-body simulation to demonstrate that the existence 
of the Sloan Great Wall and a void complex in the SDSS region 
is perfectly consistent with the ACDM model. Considering all 
these results, it should be prudent for one to claim that the recent 
discoveries actually contradic t the co smological principle of ho¬ 
mogeneity. In fact, iNadathuJ (1201 3h suggested that the homo¬ 
geneity scale is an average property. It is not necessarily affected 
by the discovery of a single large structure. It implies that very 
massive structures in the relatively local Universe could conceiv¬ 
ably be present. 

Despite the scattering results, objects like quasars and LQGs 
and events like GRBs provide new tools to study the matter dis¬ 
tribution and any possible power excess in the intermediate- and 
high-redshift Universe. We hope that the next generation of sky 
surveys will offer excellent prospects for clearing up the per¬ 
plexities between the observations of large-scale structures and 
the standard excursion set theory of structure formation. 
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